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Abstract 



The linear and cubic dynamic susceptibilities of solid dispersions of nano- 
sized maghemite particles have been measured for three samples with a volume 
concentration of magnetic particles ranging from 0.3 % to 17 %, in order to 
study the effect of dipole-dipole interactions. Significant differences between 
the dynamic response of the three samples are observed. The dynamic suscep- 
tibilities of the most dilute sample compare reasonably well with an existing 
theory for the linear and cubic dynamic susceptibilities of an assembly of 
noninteracting, uniaxial magnetic particles. The nonlinear dynamic response 
of the most concentrated sample exhibits at low temperatures some of the 
features observed in a Ag(ll at.% Mn) spin glass. 

I. INTRODUCTION 

The magnetic response of ensembles of single-domain magnetic particles has been in focus 
since the pioneering work of Neel and Brown.0'1 In the case of noninteracting particle systems, 
much theoretical wor have been devoted to the linear response. Experimental results 
from magnetic relaxation and ac susceptibility measurements!'! reasonably well support the 
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now existing models. 

The linear response of interacting magnetic particle systems is more involved. It has re- 
cently been shown that dipole-dipole interactions introduce collective behavior, as evidenced 
by the appearance of magnetic aging and a significantly broadened magnetic relaxation.i'i 
Moreover, the dynamics of a magnetic particle system of monodispersive nature is indicative 
of critical slowing down at a finite temperature,0 implying the existence of a low temperature 
spin-glass-like phase. 

Much less work have focused on the nonlinear response of magnetic particle systems. 
Bitoh and co-workers0ll! studied experimentally the third harmonic of the ac magnetization 
to obtain information related to the cubic dynamic susceptibility. However, the ac field used 
in this study was comparably large,El which makes it possible that the results for the third 
harmonic of the magnetization are contaminated with higher order susceptibility terms. In 
order to obtain theoretical expressions for the dynamic nonlinear susceptibility for systems 
without magnetic anisotropy, one can directly translate the expressions obtained for the 
dielectric relaxation.0 The inclusion of magnetic anisotropy is more involved and it is only 
recently that this has been done for systems with the simplest uniaxial anisotropy.0lilELl Of 
special interest is the study of the cubic dynamic susceptibility by Raikher and Stepanovjl! 
who used the Fokker-Planck equation to obtain numerically exact results on the linear and 
cubic dynamic susceptibilities. They also suggested approximate analytical expressions for 
these quantities. 

In the case of interacting magnetic particle systems, experimental work related to the 
nonlinear magnetic response are even more scarce.0'0 However, in a recent study it was 
shown that the equilibrium nonlinear response of an interacting magnetic particle system of 
monodispersive nature indicates that the cubic equilibrium susceptibility diverges at a finite 
temperature,^ thus providing further evidence for a low temperature spin-glass-like phase. 
To the best of our knowledge, no experimental work has been reported for the dynamic 
higher order susceptibility terms of interacting particle systems. 

In this paper, we report a study of the linear and cubic dynamic susceptibilities of 



a magnetic particle system consisting of nanosized maghemite particles. The effects of 
interparticle interactions are investigated by studying three samples with different volume 
concentrations of magnetic particles. The results for the most dilute sample are compared 
with the expressions for the susceptibility proposed by Raikher and Stepanov. The dynamic 
response of the most concentrated sample is compared with the dynamic response of a 
Ag(llat.%Mn) spin glass. 



II. EXPERIMENTAL 



The experiments were performed on three samples consisting of nanosized maghemite 
(7-Fe2C>3) particles, with a typical particle diameter of 6 nmEl and almost spherical par- 
ticle shape (observed in TEM analysis) .0 The particles were suspended in a hydrocarbon 
oil and coated with a surfactant layer preventing the particles from agglomerating. Since 
the measurements were performed at low temperatures, the oil was frozen and the par- 
ticles fixed randomly in the sample. Three samples with different volume concentration 
of particles, ~0.3%, 3%, and 17%, were used in order to study the effects of magnetic 
dipole-dipole interactions. The lowest concentration should be sufficiently low to make the 
dipole-dipole interactions between the particles negligible. Experiments were also performed 
on a Ag(llat.%)Mn spin glass sample exhibiting long range spin-spin interactions of RKKY 
type. . The spin glass sample was prepared by melting pure Ag and Mn together at 1000 °C 
in an evacuated and sealed silica tube. After annealing the sample at 850 °C for 72h it was 
quenched in water to room temperature. 

Two different experimental equipments have been used: (i) A commercial ac- 
susceptometerlH was used to measure the 1st and 3rd harmonics of the magnetization for 
different ac-field amplitudes in the range 100 - 2000 A/m. (ii) A non-commercial low-field 
superconducting quantum interference device (SQUID) magnetometer^ was used to perform 
sensitive studies of the frequency dependence of the linear and cubic dynamic susceptibilities. 
Frequencies in range of 2 Hz to 200 Hz were used. 
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The magnetization, M, can be expanded with respect to an applied field H as 

M = X H + X3H 3 + X5 H 5 + ---, (2.1) 

where x is the linear susceptibility and X3 the cubic susceptibility. The dynamic susceptibil- 
ity can be probed by applying an ac field, H = ho cos (ut). The linear susceptibility is then 
obtained from the magnetization measured at the fundamental frequency as x( u ) = M^/ho 
and the cubic dynamic susceptibility is obtained from the third harmonic of the magnetiza- 
tion as jXsiu) — M^/hQ. These equations are only valid if the applied ac field is sufficiently 
small, so that contributions from higher order susceptibility terms in the expansion of M w 
and are negligible. For the measurements of the linear and cubic dynamic susceptibil- 
ities we used different ac-fields to ascertain that no mixing with higher order susceptibility 
terms occured. 



III. THEORETICAL BACKGROUND 

A. Equilibrium susceptibilities 

The potential energy of a single-domain particle in an external magnetic field is the sum 
of the Zeeman energy and the anisotropy energy. For a particle with uniaxial anisotropy the 
total magnetic potential is given by 

U = -D(e-h) 2 - n M s V(e- H), (3.1) 

where D is the anisotropy energy barrier and H is the external filed, while n and e are the 
unit vectors along the anisotropy axis and magnetic moment, respectively. In the case of 
volume anisotropy D = KV, where V is the particle volume and K is the volume anisotropy 
constant, and for surface anisotropy D = K$S, where S is the particle surface and K$ is 
the surface anisotropy constant. 

The linear and cubic equilibrium susceptibilities can be derived from the partition func- 
tion associated with the magnetic potential. For a mono-dispersive system with randomly 
distributed anisotropy axes, the linear and cubic equilibrium susceptibilities are given bySil 
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respectively, where 

S 2 (ct) = Z - x ((7) y P 2 (z) exp(az 2 )dz. (3.3) 

Here i^^) = |(3-2 2 — 1) is the second Legendre polynomial, z = (e • n), a = D/k B T, and 
Z = j\ exp(az 2 )dz is the partition function in the absence of an applied magnetic field. 

B. Dynamic susceptibilities 

Recently, Raikher and Stepanovll! studied theoretically the linear and cubic dynamic sus- 
ceptibilities of noninteracting ensembles of single-domain particles with uniaxial anisotropy, 
solving numerically the Fokker-Planck equation in the overdamped case. It is difficult to 
derive analytical expressions for these quantities, but for a particle system with randomly 
distributed anisotropy axes, they suggested the following simple formulae for the linear 
dynamic susceptibility, 



X( T , u) = x 



cq 



(3.4) 



_3(l + i0JT) 3 

This expression has been shown to be a good approximation to the exact linear dynamic 
susceptibility for frequencies below the ferromagnetic resonance frequency regime (see, for 
instance, Ref. |23|). In the above expression r is the relaxation time, for which various 
analytical expressions have been suggested (see, for instance, Ref. |4|). In the moderate to 
high-energy barrier case, the relaxation time is approximately given by an Arrhenius law 
as r = t exp(D//c B T), where r is a constant. Raikher and Stepanov also proposed an 
expression for the cubic dynamic susceptibility 

(rr \ eg (1 ~ /„ ,N 

(1 + IUT)(1 + dlUJT) 

which they assert to be a good approximation of the numerically exact result in the low 
frequency regime. From the above equations we can see that for lot 1, e.g., for high tem- 
peratures, the dynamic susceptibilities reduce to the equilibrium susceptibilities Eq. 



IV. RESULTS AND DISCUSSION 



The linear susceptibilities for the three maghemite samples measured at the frequency 
lv/2-k = 125 Hz and with an ac-field amplitude of 100 A/m are shown in Fig. |l|. Dipole-dipole 
interactions shift the susceptibility peaks to higher temperatures, lower the magnitude of the 
peaks, but produce a higher equilibrium susceptibility. Figure |2] shows the corresponding 
cubic susceptibilities measured with an ac-field amplitude of 400 A/m. Similar to the case 
of the linear susceptibility, dipole-dipole interactions shift the susceptibility peaks to higher 
temperatures and reduce their magnitudes. The sample with the highest concentration of 
particles has a second positive peak at low temperatures in both the real and imaginary 
components of the cubic susceptibility. 

The results for the most dilute sample and the most concentrated sample will be dis- 
cussed separately below. For the most dilute sample, the linear and cubic susceptibilities 
are compared to the theoretical expressions for noninteracting particle systems discussed in 
Sec |T|. The sample with the highest particle concentration has been reported to show col- 
lective behavior at low temperatures J Its cubic dynamic susceptibility is therefore compared 
with xs obtained for the spin glass sample. 



A. Noninteracting particles 

To compare the measured linear and cubic susceptibilities with theoretical expressions, 
the polydispersivity of the particle system needs to be taken into account. The conventional 
approach is to choose a trial volume distribution and determine its parameters by fitting 
theoretical curves to experimental data. The log-normal distribution and the gamma distri- 
bution have been shown to work well for a variety of nanoparticle systems. The log-normal 
distribution is given by 

s00 = -i «p(-!=2MV (4.i) 



TkVs \ 2s 2 
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where Vq is the median particle volume and s is the logarithmic standard deviation. The 
gamma distribution is given by 

9{V) = TV-, 1 oM, (tf) exp(-VyVo), (4.2) 



V(1+P)V Wo, 

where T(-) is the gamma function and Vq and f3 are parameters related to the mean particle 
volume and standard deviation. 

The susceptibility of a polydispersive system is calculated as Xpoiy = Jy XmonogvdV , 
where gv is the volume weighted volume distribution and Xmono is either the linear or the 
cubic susceptibility given by Eq. (|3.4|) and Eq. Q3.5|) , respectively. In the theoretical model 
the anisotropy is assumed to have uniaxial symmetry. In other work on 7-Fe203 nanoparti- 
cles, it has been shown that the dominating magnetic anisotropy is uniaxial and of surface 
type.0^1 In these studies particles exhibiting almost spherical particle shape and with di- 
ameters in the range 4.8 - 10 nm were investigated]^ which correspond well to the particles 
investigated in the present study. 

The temperature dependence of the saturation magnetization at low temperatures has 
earlier been modeled as M S (T) = M S (0)(1 - 1.8 x 1(T 5 T 3 / 2 ), with M s (0) = 4.2 x 10 5 A/m.0 
The Nelder-Mead simplex method was used to perform nonlinear fitting to the measured 
linear and cubic susceptibilities for frequencies in the 2 Hz to 200 Hz range. The fitting 
parameters were: the anisotropy constant, the preexponential factor in the expression for the 
relaxation time (r ), and the two parameters from the volume distribution. The parameters 
obtained from simultaneous fits to \ an d X3 are presented in Table |. (We do not present the 
combination surface anisotropy and volume distribution modeled by a gamma distribution, 
since it gave unphysical fitting parameters.) The anisotropy constant in the case of surface 
anisotropy K$ = 2.3 x 10~ 5 Jm -2 is similar to the value K$ = 2.7 x 10~ 5 Jm -2 obtained in 



Ref. pTL The different volume distributions give a typical particle diameter of 6 ± 1 nm in 



good agreement with the typical particle size of 6 nm observed in TEM studies. @B The 
results of the simultaneous fits to x an d X3 are shown in Fig. for the linear susceptibility 
and in Fig. |4| for the cubic susceptibility. We can see that the best fit is obtained for a volume 



distribution modeled by a gamma distribution. However, the differences between the three 
cases are not sufficiently large to make any conclusions about the origin or the uniaxial 
anisotropy or the functional form of the volume distribution in this sample. We have also 
performed the fitting with the expression for the relaxation time proposed by Cregg et al. 
in Ref. but the change of the quality of the fit was negligible. 

The discrepancies between the calculated and the measured dynamic cubic susceptibili- 
ties may have several possible origins. One obvious origin is that the models of the polydis- 
persivity of the system used here are too simple and that the real situation is more involved. 
It may be that the symmetry of the magnetic anisotropy in the experimental system is 
different from uniaxial. Such a difference need not originate from the magnetocrystalline 
anisotropy but can also be due to the geometrical shape of the individual particles. The 
lack of theoretical work where a magnetic anisotropy different from the uniaxial case is as- 
sumed makes it difficult to predict what effect the symmetry of the anisotropy will have 
on the nonlinear dynamic susceptibilities. Another possible cause for the observed differ- 
ences is linked to the details of the Fokker-Planck equation used in the theoretical work of 
Raikher and Stepanov,0 who studied the overdamped case (the damping constant, in the 
dynamical Landau Lifshitz equation from which the Fokker-Planck equation is derived, is 
assumed to fulfill a ^> 1). However, a more realistic value of the damping constant for 
magnetic nano-particles is a < 1.0 The form of the dynamic linear susceptibility curves 
is not very sensitive to the value of the damping constant (it essentially changes the value 
of To), but the same does not necessarily hold for the nonlinear response and a relatively 
more pronounced dependence of Xz{ u ) on the damping parameter cannot, in principle, be 
ruled out. Taking all these considerations into account, we can conclude that measurements 
of both the linear and cubic susceptibilities potentially allow for extracting more detailed 
information concerning the intrinsic properties of real particle samples. 
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B. Interacting particles 



Figures |l] and |2| clearly show that dipole-dipole interactions have a strong influence on 
the magnetic relaxation of the particle system. Some recent experimental studies show that 
strongly interacting particle systems approach a spin glass phase at low temperatures. For 
example, spin glass characteristics such as: aging effects,^ critical slowing down,0 and a 
divergence of the cubic equilibrium susceptibility^ have been reported in such systems. 

In a model often used to analyze the dynamics of interacting particle systems,^ the effect 
of interparticle interactions is accounted for by shifting the energy barrier distribution to 
higher energies and thereby increasing the relaxation times. Some features shown in Figs. [I] 
and ^ such as a shift of the maxima in x(T, uj) to higher temperatures and the appearance of 
weak positive peaks in X3(T, uj) at low temperatures, can in principle be explained by a signif- 
icant increase of the individual particle relaxation times. However, for the here investigated 
particle system, it has previously been shown that the low temperature magnetic relaxation 
for the most interacting sample resembles that observed for archetypal spin glasses.! At tem- 
peratures below 45 K, the relaxation time spectrum broadens significantly and the magnetic 
relaxation is significantly different from that of the noninteracting sample. This together 
with the observed aging effect! give strong indications that the low temperature magnetic 
relaxation is dominated by collective particle dynamics. With increasing temperature, how- 
ever, the time scale of the collective dynamics is gradually shifted to shorter time scales and 
the slow magnetic relaxation remaining is due to single-particle relaxation of a comparably 
small number of large particles. Still, the relaxation of these large particles is influenced 
by small particles surrounding them - relaxing particles will experience a magnetic field 
from neighboring polarized superparamagnetic particles, increasing the energy barriers and 
thereby the relaxation times. This gradual change of the magnetic relaxation with increas- 
ing temperature, from collective to single particle dynamics, obstructs the observation of 
critical slowing down in this particle system. For the same reason, it will not be possible 
to experimentally observe a divergent cubic equilibrium susceptibility. Nevertheless, the 
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observation of collective behavior of the magnetic relaxation at low temperature makes it 
interesting to compare the behavior of the most interacting sample to that of a typical spin 
glass with long range interactions. 

Figure |5] shows the susceptibility of a Ag(ll at.% Mn) spin glass measured at the fre- 
quency u/2tt = 125 Hz and with an ac field amplitude of 1600 A/m. The real part of the 
linear susceptibility has a sharp cusp at about the same temperature at which there is a 
sharp rise from zero of the imaginary part. The two components of the cubic susceptibility 
have sharp negative peaks at high temperatures followed by broad positive peaks at low 
temperatures. Before comparing this behavior with that of the most concentrated parti- 
cle system, we will point out some fundamental differences between an interacting particle 
system and a spin glass. 

In a particle system the single particle relaxation time has an exponential temperature 
dependence and in addition depends upon the particle volume and the anisotropy constant, 
while in a spin glass the individual spins have a relaxation time in the order of 10" 12 s to 
10~ 14 s. In a spin glass the slow dynamics is only due to collective phenomena and thus the 
imaginary component of the susceptibility disappears when the maximum relaxation time 
of the spin system becomes shorter than the observation time 1/u. The dynamics of the 
investigated particle system is, as discussed above, dominated by collective particle behavior 
at low temperatures, while at higher temperatures the slow dynamics are dominated by single 
particle relaxation. When comparing the two systems, we are thus restricted to the low 
temperature regime where both systems display collective dynamics. For the investigated 
interacting particle system, this implies temperatures T < 45 K.i 

The linear and cubic susceptibilities of the most concentrated particle sample (see Figs. 
[I] and 0) display features much broader in temperature as compared to the corresponding 
features observed for the spin glass sample. This can be explained as a combination of 
two effects: the polydispersivity of the particle system and that the maxima in the linear 
susceptibility as well as the negative peaks in the cubic susceptibility are partly due to single 
particle relaxation. It is tempting though to attribute the low temperature positive peaks in 
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the interacting particle sample to collective dynamics, since also the spin glass shows broad 
positive peaks in the cubic dynamic susceptibility at low temperatures. This feature is found 
at temperatures T < 40 K, i.e., in the temperature range where the aforementioned study 
of the magnetic relaxation of the most interacting sample revealed evidence for collective 
particle behavior. 

V. CONCLUSIONS 

We have studied the linear and cubic dynamic susceptibilities of solid dispersions of 
maghemite nanoparticles with various strengths of the inter-particle interactions. We have 
found that the expressions for the dynamical susceptibilities proposed by Raikher and 
Stepanov, describe the experimental results for the most dilute sample with a reasonable de- 
gree of accuracy. Nevertheless, further developments of the theoretical modeling (including 
other symmetries of the anisotropy energy, effects of the damping parameter, etc.) would 
be desirable to reduce the gap between theory and experiments. 

Concerning the features observed in the most concentrated sample an explanation in 
terms of single particle dynamics with equilibrium susceptibilities and energy barriers modi- 
fied by the interactions cannot be excluded only on the basis of the present study. However, 
previous evidence of collective magnetic behavior in the same system strongly suggests this 
as a more consistent interpretation of the obtained results. This is also supported by the 
comparison of the results obtained for the most concentrated particle system with the non- 
linear response of an archetypal spin glass. 
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TABLE I. Parameters describing the particle system obtained by simultaneous fitting to the 
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Volume distribution Vq [nm 3 ] Anisotropy constant To [s] 

Log-normal 230 s = 0.68 K = 1.4 x 10 4 Jm~ 3 

Log-normal 200 s = 0.78 K s = 2.3 x 10~ 5 Jm~ 2 

Gamma 210 (3 = 0.34 K = 1.2 x 10 4 Jm~ 3 
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FIG. 1. The real (a) and the imaginary (b) parts of the linear susceptibility vs. temperature, 
for the three samples with different volume concentration of magnetic particles. lo/2it = 125 Hz 
and ho = 100 A/m. 
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FIG. 2. The real (a) and the imaginary (b) parts of the cubic susceptibility vs. temperature 
measured at the angular frequency 3u>, with ui/2tt = 125 Hz, for the three samples with different 
volume concentration of magnetic particles, ho = 200 A/m. 
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FIG. 3. The linear susceptibility vs. temperature for the sample with a volume concentration 

of 0.3%. Solid lines represent the calculated susceptibility using volume anisotropy and the volume 
distribution taken as a gamma distribution, dashed lines show the corresponding result but with 
the volume distribution taken as a log-normal distribution and dashed-dotted lines represent the 
calculated susceptibility using surface anisotropy and the volume distribution taken as a log-normal 
distribution. 
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FIG. 4. The cubic susceptibility vs. temperature for the sample with a volume concentration 

of 0.3%. Solid lines represent the calculated susceptibility using volume anisotropy and the volume 
distribution taken as a gamma distribution, dashed lines show the corresponding result but with 
the volume distribution taken as a log-normal distribution and dashed-dotted lines represent the 
calculated susceptibility using surface anisotropy and the volume distribution taken as a log-normal 
distribution. 
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FIG. 5. The real (a) and the imaginary (b) parts of the cubic dynamic susceptibility vs. tem- 
perature measured at the angular frequency 3uj for the Ag(ll at.% Mn) spin glass sample. The 
linear susceptibility measured at oj/2tt = 125 Hz is shown in the insets, ho = 1600 A/m. 
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